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Let J 2 be the set of n X n complex matrices A = (aij) such that aj^a^ . . . aj r j 1 = for all r such 
that 3 < r ^ n and all distinct j u ji, . . ., jr. Then many properties of this set are given, which may be 
regarded as generalizations of the properties of the set of triple diagonal matrices. 
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1. Introduction 

In a previous paper, [l], 1 we introduced the class 
of matrices ^ 2 as a generalization of the class J' of 
Jacobi matrices. We also developed a few of the 
spectral properties of certain types of matrices in 
^2. In [2] we have shown why the members of J^ 2 enjoy 
many special properties not shared by matrices with 
more complex structure. 

Several years prior to the work done in [1] and [2] 
S. Parter in [3] obtained a spectral theorem for certain 
elements of ^ 2 by an ingenious argument based upon 
the theory of linear graphs. 

Finally, in a recent paper [10], Quirk and Ruppert 
obtained very deep stability results for a subset of 
^2 quite different from that considered by Parter. 
The theorems in [1] overlap with the work of Quirk 
and Ruppert as well as with the work of Parter. 

The various results of the papers [1], [2], [3], and 
[10] point to the desirability of constructing a compre- 
hensive theory of the properties of elements of ^ 2 . 
The purpose of the present paper is to initiate such a 
study. We will show that a remarkable number of the 
properties of Jacobi matrices can be generalized to 
the elements of J>%. 

It turns out that the class J r 2 contains two subclasses 
of particular importance. We shall refer to the elements 
of these classes as matrices of semi-positive type and 
semi-negative type respectively. The semi-positive 
matrices occur in a variety of problems in the physical 
sciences, particularly in the small vibrations of 
mechanical systems. The semi-negative matrices 
occur in the theory of qualitative stability and find 
their primary application in the social sciences, 
particularly in economics. 

With these applications in mind we have concen- 
trated mainly upon the semi-positive and semi-negative 
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elements of £% although these two classes certainly 
do not exhaust the possibilities. 

The paper is organized into several parts as follows: 
In section 2 we set forth the basic definitions and 
develop what may be called the combinational proper- 
ties of the elements of ^ 2 . This section is not restricted 
by any considerations of positivity or negativity. 
The same is true of section 3 in which we give a variety 
of determinant formulas. Some of these results gen- 
eralize known formulas for Jacobi matrices, but the 
problem for elements of ^ 2 is more difficult and some 
problems remain open. In particular, our formulas for 
the exterior p-th power of elements of ^ 2 are not suffi- 
ciently well worked out to enable us to decide which 
matrices have exterior p-th powers of positive cyclic 
type. On the other hand many of these formulas hold 
for quite general matrices. 

Section 4 consists of a study of semi-positive ele- 
ments of ^ 2 . Here we concentrate mainly upon the 
interlacing and spectral multiplicity, 
special, aspects of the semi-positive 
more properly reserved for a general 
treatment of certain mechanical problems where such 
questions naturally arise. We intend to treat these 
matters shortly in another paper. 

Finally, section 5 consists of a few miscellaneous 
properties of semi-negative matrices in ^ 2 . The prob- 
lems here are more difficult than those associated 
with the semi-positive matrices because they center 
around the problem of stability and this in turn seems 
to be an inherently difficult problem. We have been 
unable to give an essentially simpler proof of the 
fundamental stability theorem of Quirk and Ruppert 
[10]. This being so, we have confined ourselves to 
stating the theorem and to deriving those few results 
which follow readily from the general determinant 
formulas of section 3. 

A basic problem remaining unsolved is that of deter- 
mining for the elements of ^ 2 how the eigenvalues of 
A, or at least the real parts of the eigenvalues, vary 
with the elements along the principal diagonal of A. 
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A solution to this problem in combination with lemma 
5.2 might yield an elegant proof of half of the Quirk- 
Ruppert theorem. Another problem not treated here 
is the structure of the eigenvectors for the semi- 
positive and semi-negative matrices. For certain im- 
portant subclasses of semi-positive matrices this 
question is most strongly motivated by problems in 
mechanics and we intend to deal with it later. Parter 
[3] has some results about eigenvectors, but a glance 
at the extensive results available for the Jacobi 
matrices makes it clear that much more extensive 
theorems should be available for J> 2 . As of the present 
writing nothing is known about the eigenvectors of 
semi-negative elements of J> 2 . 



2. Notation and Fundamental Preliminaries 

We deal with n X n matrices over the real or complex 
field, n > 2. A = (a,*)? is a Jacobi matrix, A e J>, if 
ajk = whenever \j — k\ > 1. To define the class J> 2 
we recall that an r-cycle in A is an object of the form 



a(ji 



• h) = aj 



(2.1) 



where the indices in J = (/1 • ■ • jr) are distinct ele- 
ments of the set S n — (1 . . . n). The elements djj on 
the principal diagonal of A are the 1-cycles of A. 

Definition 2.1: The matrix A = (a jk y{ is an element 
of J> 2 if every r-cycle of A of length r > 2 is zero. 

Definition 2.1 leads at once to a formula of funda- 
mental importance in any study of the properties of 
the elements of J> 2 . Let J = (j\ . . . 7V) be any set of 
distinct indices from S«; then we denote by A[J] (or 
A\}\ • • • jr]) the principal submatrix of A in rows and 
columns J. We use both the symbol d to denote the 
determinant function and the Bourbaki notation. 
Thus, 

d{A[J])=Aj = d(A\j, . . . j r ]). (2.2) 

We also write A(k) for A - kl, Aj(K) = d(A[J; A]). For 
each r=2, 3, . . ., n, each sequence J= (j\ . . . j r ) 
of distinct indices, and each 1 =^ p ^ r, we have 



Aj(k)=(a(jp) ~ k)Aj- jp (k) - £ UjtiMj - (jj,) (X). (2.3) 

q=l 
q^p 

In this formula and subsequently we denote by J — H 
(or J — (h\ . . . h s )) the set of indices remaining in J 
after deleting the indices in H. 

In [1] we defined the class ^ 2 directly by a formula 
having the form (2.3). We were led naturally to it by 
the corresponding formula connecting the principal 
minors of Jacobi matrices. 

The key to our method consists in first showing that 
practically every determinental formula normally 
associated with the Matrix A can be expressed in 
terms of determinents of principal submatrices and 
certain special products of elements of A which are 
easily studied. The new formulas obtained in this way 
are then systematically exploited. We have already 



done this with the basic determinant formula intro- 
duced in [2] and the results obtained there are freely 
used below. Formula (2.3) is at the same time a special 
case of the classical formula of Cauchy and a special 
case of the general formula of [2] for an arbitrary square 
matrix. 

To expedite the study of various products in a matrix 
A we introduce the following concepts. 

Definition 2.2: The product 



a(J) = a(ji • • ■ Jp) = a J 



Jp-iJp' (2-4) 



(/1 . . . jp) distinct, is called a chain of length p — 1 
in A. 

In dealing with chains we shall often refer to 
a(aj(3) = a(aji . . . j p f$) as a chain from a to /3. 
The number of distinct chains in A of length p — 1 is 



pll I. A complete theory of chains and cycles of 

V 

A is not essential for the study of J> 2 so we limit 
ourselves here to the few basic concepts and ele- 
mentary results needed. 

Definition 2.3: The product 



aC/)=aC/)a 



l j\jp 



(2.5) 



is called the closure of a(J). We also refer to a(J) 
as a closed chain of length p in A. The cycle 

«l/)=«Vi 

will be referred to as the completion of a(J). 

Note that a chain of length 1 is closed. All other 
closed chains have length ^ 3. In the sequel when we 
refer to closed chains we shall exclude the closed 
chains of length 1. 

Clearly there is a unique closure and a unique com- 
pletion to each chain in A. Conversely each closed 
chain is the closure of a unique chain so that there is a 
1-1 correspondence between chains of length p — 1 
and closed chains of length p. On the other hand a 
cycle of length p is the completion of each of the p 
chains of length p — 1 obtained by striking out any 
term in a (J). 

A chain a(J) +0 will be called maximal in A if both 
its closure and its completion are zero. Clearly we 
may also characterize the class J 2 by the statement 
that for A € J> 2 no chain of length greater than one has a 
nonzero completion in A. Recalling (see [1], [2], or 
[3]) that A is combinatorially symmetric if a jA - 4= 
implies a^-j + O, we see that the existence of nonzero 
closed chains in A e J> 2 is a measure of the extent to 
which A fails to be combinatorially symmetric. On 
the other hand we shall see that A € J> 2 is irreducible 
(indecomposable) only in the presence of combinatorial 
symmetry. Hence the existence of one nonzero closed 
chain in A e J^ 2 implies A is reducible. This furnishes 
us with an occasionally useful criterion for recognizing 
reducible elements of ^ 2 . 

Definition 2.4: The set Q consists of all A e J 2 
which are also irreducible (indecomposable). 
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In [1] we called the elements of Q quasi-Jacobi 
matrices. The paper [3] of Parter deals with a proper 
subset of Q. Some of our results will be valid only for 
elements of Q; others will be valid for all elements 

Of/2- 

Theorem 2.1: Let A e J 2 . A e Q if and only if A is 
combinatorially symmetric and has exactly n — 1 
nonzero elements above the principal diagonal. 

PROOF: We confine ourselves here to proving the if 
portion of the theorem since the only if portion was 
proved in [1] (see also [10]). Suppose, for contradiction, 
that the set S n can be divided into disjoint subsets 
/ = (/i • • • jp)i J' = (j)>+i • • • jn) such that every ele- 
ment having the form ajk with j ej, he J' is zero. By 
combinatorial symmetry the elements a^j are also 
zero. The principal submatrices A[J] and A[J f ] are also 
combinatorially symmetric. Either A[J] contains more 
than p — 1 nonzero elements above the principal diago- 
nal or A[J'] contains more than n—p — l nonzero 
elements above the principal diagonal. For in the 
contrary case A can have no more than n — 2 nonzero 
elements above the principal diagonal. By theorem 1 
of [1] it follows that A >k J>t yielding the desired contra- 
diction. Thus the conditions imply A is irreducible. 

We may observe now on the basis of this theorem 
and the previous remarks that A € Q if and only if 
every chain in A of length ^ 2 is maximal. 

For many purposes it is desirable to know the rela- 
tive positions of the elements of A € J> 2 which are 
zero and the nonzero elements of A. Of course, each 
nonzero off diagonal element of A is a chain of length 
one. The existence of nonzero chains of greater length 
in A does imply information on the location of zeroes. 
The following result is almost obvious. 

Lemma 2.2: Let a(ii . . . i p ) be a nonzero chain in 
AeJf-2 of length p^2. Then every element a„# of A 
such that a = i q , (3 — i r with l^r^q— l=^p is zero. 

The lemma implies that each nonzero chain of length 
p — > 2 in A e J>2 forces (p — l)(p — 2)/2 elements of 
A to be zero. If A e Q the number of zeroes becomes 
(p— l)(p — 2). It is instructive to look at the two 
extreme cases below (the x's denote the location of 
nonzero elements), 
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A u A 2 eQ. A x has the chain a( 124653) 4= of length 
5. By lemma 2.2 this nonzero chain locates all of the 
zero elements of A\. On the other hand A 2 clearly has 
no nonzero chain of length greater than 2. 

Here is a result which is complimentary to lemma 2.2 
and which we again give without proof. 

Lemma 2.3: Let A e Q and suppose that for fixed 
1 ^ i ^ n the elements a ih , . . ., a,j , ji 4= i, . . ., 
j q =t= i, are all different from zero. Then the 

(q-l)(q-2)/2 

elements a jpjr with l^r^p— l=^q of A are zero, 
and the symmetrically placed elements are zero. 

These simple ideas may be put to work to further 
classify the elements of Jf 2 with aid of 
Definition 2.5: The matrix AeQ belongs to the set 
Q n ~ l if there exists a nonzero chain of length n — \ 
in A. 

Theorem 2.4: The matrix A e Q n_1 if and only if 
there exists a sequence A°=l, A 1 , . . ., A n =A of 
principal submatrices of A with A J a matrix of order 
j, j = 0, . . ., n, and A j_1 a principal submatrix of 
A j , satisfying the recurrence formula 

d(A cl )=a(j q )d(A <l - 1 )-a(j q j a - 1 )d(A tl -2), (2.7) 

q = 2, . . ., n, where the products a(j q j q -i) 4= 0. 
PROOF: Suppose first that A eQ" ' and let 



«0'i 



• • Jn) 



(A , 

«0\ • 
a'(j\ . 
n-l 
cipal 



be a nonzero chain of length // — 1 in A. Set A°=\, 
Jq = Ui • • • jq] f° r each q=l, . . ., n, and set 
A Q =A[J Q ]. Expand d{A q ) relative to the gth column. 
By virtue of lemma 2.2 the elements aj , = for 
r ^ 2. On the other hand, we are given that aj j =1= 0. 
By combinatorial symmetry the only nonzero off- 
diagonal element in the qth row of A q is a,, . For- 
mula (2.7) follows immediately for q ^ 2. 

To prove the converse, suppose the sequence 

. . ., A 11 ) exists satisfying (2.7). The chain 

• -> jn) + and also the transposed chain 

. .j n ) = a(j n . . . j\) 4=0. Thus A has at least 

nonzero off-diagonal elements above the prin- 

diagonal and at least n — 1 below. Since AeQ, 

AeQ" 1 . 

As an application of this result consider the matrix 
A above. A sequence of principal submatrices satis- 
fying the formula (2.7) is the following: 

A 1 = A[l],A* = A[12] 9 A* = A[m'], 

A 4 = A[1246],A* = A [12456], A« = A [123456]. 

The sequence A 1 =A[3], A 2 = A[35], A 3 = A[356], 
^ 4 = ^[3456], A* = A [23456], A« = A also satisfies 
such a recurrence formula. It is clear that in general 
there are exactly two sequences in AeQ nl satisfying 
the formula (2.7), one sequence running in the oppo- 
site direction from the other. 
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Theorem 2.4 can be put a little differently by saying 
that AeQ" 1 if and only if there is a nested sequence 
of principal minors of A satisfying a second order dif- 
ference equation. 

Here is an alternative important characterization 
of Q n ~\ 

Theorem 2.5: Let AeQ. Then AeQ nl if and only 
if no row (or column) of A contains more than two non- 
zero off-diagonal elements. 

PROOF: Assume first that AeQ^ lX but that row i 
contains the nonzero elements a^, a\ q , a* r , 

(p =1= i , q =t= i , r =\= i ) . 

A contains the nonzero chains a(i\ . . . i n ) and 
a' (ii . . . i n ). The product 

a(U . . . i n )a'{U • • . in) =N 

contains all of the nonzero off-diagonal entries of A. 
Each subscript appears at most twice in a(U . . . i n ) , 
once as a row subscript and once as a column sub- 
script, and the same is true of a' (i\ . . . i n ). But if 
dip, aiq, air are nonzero so are a P i, a q i, a r i and the 
product a(i\ . . . i n ) must contain all six of these 
elements so that the index i appears six times in the 
product. This contradiction establishes the only if 
portion of the theorem. 

Assume next that AeQ and the condition is satis- 
fied. Denote by a(U • . • inin • • • ii ) =♦= the prod- 
uct of the /i-l nonzero 2-cycles of A. Each index 1 
through n appears in this product at most twice as a 
row subscript and at most twice as a column sub- 
script, and some pair of indices appear just twice. 
This follows from the irreducibility of A. That AeQ" 1 
is now clear. 

Let us now complete our remarks on Q' 1 ' 1 by the 
following result: 

THEOREM 2.6: Let AeQ' 1 " 1 , then there exists a permu- 
tation matrix P such that the similarity PAP trans- 
forms A into a Jacobi matrix. 

PROOF: We have merely to observe that the permu- 
tation such that T(j Q ) = q will convert the recurrence 
relation (2.7) into the usual relation for a Jacobi 
matrix. 

For the example A, above the permutation 

o-(123456) = (356421) 

will do the job. 

Theorem 2.6 shows that the elements of Q n ~ x are 
essentially Jacobi matrices. We shall exploit this 
fact below. 



ties of the elements of Q. Therefore we shall briefly 
remark on these connections first. 

Definition 3.1.: Let A= (ajk)i be combinatorially 
symmetric. By the graph of A, G(A), we mean a set 
{pi, . . ., pn] together with certain distinguished 
sets of pairs {pt, pj} corresponding to the elements 
aij =f= of A. The points pi, . . . , p n are called the 
vertices of the graph and the pairs {p*, pj] the arcs of 
the graph. 

Note that this graph G(A) is only defined here for 
combinatorially symmetric matrices, hence it is not 
the directed linear graph used, for example, by Varga 
in [13]. 

G(A) is called a tree if it has Betti number zero 
and is connected. The next result is a direct conse- 
quence of the concepts of section 2. 

Theorem 3.1: Each of the following two conditions 
is both necessary and sufficient in order that G(A) be 
a tree. 

(1) AeQ, 

(2) For each pair of indices i =t= j there is exactly one 
nonzero chain a(i J k), J = (j, . . ., j p ). 

The condition (2) of the theorem turns out to play 
a crucial role in our examination of the exterior p-th 
power of A. 

Clearly there is a 1-1 correspondence between 
trees and the elements of Q. Parter exploits this fact 
systematically in obtaining his results. 

We turn next to an examination of the matrices 
AM, p = 2, . . ., n— 1, A*A being the exterior p-th 
power of A. We remind the reader that if the row vec- 
tors of A are a 1 , . . ., a n , then the row vectors of 
A'A are just the exterior products 

a j u\a j zA . . . Aah 

arranged in lexicographic order. Thus A*A=A, 
A n A = d(A), and A"~ l A is, except for certain signs, 
the cofactor matrix of A. (The matrix A P A is usually 
called the p-th compound matrix of A in the older 
literature.) One may pass from the matrix A n ~ x A to the 
cofactor matrix, cof A, by multiplying they, Ar-th ele- 
ment A n ~ A A by (— iy +k . We start with cof A, since it is 
the simplest case. 

For fixed a, /3 we shall denote by A a p the algebraic 
cofactor of a a p so that cof A has the elements /!«£, 
a, 0=1, . . ., n. 

Lemma 3.2: If a 4= /3, we have 



3. Various Determinant Formulas 

In this section we collect several determinant 
formulas together with a few of their more immediate 
applications. Some of these formulas have interesting 
connections with the linear graph of A, G(A), as intro- 
duced by Parter in his study of the spectral proper- 



A a B 



a 2 « 



a/31 



Ha(3)' 



a/3n 



(3.1) 



where A( a py is the principal submatrix of A of order 
n — 2 consisting of the elements in all rows and columns 
except a and (3. 
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PROOF: By definition 



^ = (-1)^ 



a n • . 

Oa-\, 1 
#a+l, 1 

a n \ . . 



«1, £-1 



««, j8-l 



. <2i 


, 0+1 • 


• • a xn 

da- 1 
(kx+l 


. a n 


, fl+1 • 


• • dnn 



(3.2) 
a + /3 — 3 interchanges of rows and columns converts 
(3.2) into (3.1). 

We now expand the determinant (3.1) relative to 
the element a$ a using the fundamental determinant 
formula of [2]. This yields 



d a p — — [a/3ad(A( a py) 



+ 2 (-l) B " r S AjA(a\fih>y]. 

r=0 Jeilr, (a/3)' 



(3.3) 



The notation of formula (3.3) is as follows: 

(a) fi r , (a/3)' is the set of all increasing sets of r dis- 
tinct indices from S„ minus a and /3; 

(b) J' is a compliment of J relative to either one of 
the sets S n — a or S„ — /3; (thus J' has n — 1 — r 
distinct elements); 

(c) A(a\/3)(j>) is the sum over all complementary 
cycles to Aj in the matrix of A a p, (each cycle in 
the sum is an n — 1 — r cycle). 

A typical element of A (a|j8)(j'> is the sum of chains 
having the form 

a(pKa) 

where K= (/ci . . . k q ) Cj\ In particular, every chain 
in A from (3 to a appears as a coefficient in some 
term in (3.2). Thus we have 

Lemma 3.3: A a /3 = if every chain from (3 to a is 
zero. 

Lemma 3.4: Let AeQ, then for each fixed a =¥ /3 
the expansion (3.2) for A a/ s has at most one nonzero 
term. 

By lemma 3.4 formula (3.3) can be put into the form: 



Aafi=±a(l3Ja)A r 



(3.3') 



for AeQ and some suitable index set J. Therefore a 
cycle of cof A (or of A n ~ l A which has equal cycles) 
has the form: 



A\.\ 



A Vi 



= ±a{ijpi p ) 



a(iJ 1 ii)A J{ . . . Ajfr (3.4) 



Now it is clear that a(ij p i p ) . . . a{i<2j\i\) is a prod- 
uct of cycles of A on the index i\ of length ^ 2. Hence 
the product can only be different from zero if each 
cycle is of length 2. It follows that an even number 
of the chains in (3.4) must be of odd length. Referring 
back to formula (3.3) we see that the sign on the right 
side of (3.3) is therefore positive. Let us formalize 
this result. 

Lemma 3.5: Let AeQ, then every nonzero cycle of 
length ^2 of A n_1 A (or of cof A) is a product of 2- 



cycles of A and principal minors of A of the form 
(3.4) with a positive sign. 

As one might expect, the structure of the matrix 
A P A for p = 2, . . . n — 2 is more difficult to analyze. 
Here is how the formula for an arbitrary minor of A 
can be related to the principal minors of A. Let R, S 
be two increasing sets of indices of length p and sup- 
pose the number of indices in RP\S is a. Let Ak,s 
be the determinant of the submatrix of A in rows R 
and columns S and let p = p — a. After suitable in- 
terchanges of rows and columns we arrive at the 

submatrix A [/?, S] having the form, 



A[R,S] = 



a/I,*-, • 


• • °Va 


a *a*l • 


• • a »a k c 


a h k \ ■ 


■ ■ a h* a 


a Vi • 


■ ■ <Va 



a 'Vi 



ai,j. 






A[J] 



in which the index sets 

H={h x . . . h a ),K=(k, . . . k a ),J=(ji • • • jfi) 

are pairwise disjoint and A [J] is a principal submatrix. 
Clearly 

A B.B" ±A n.r 

On the other hand, repeated application of the 
determinant formula of [2] to Aii.s shows that the 
terms in the expansion of this determinant are ob- 
tained by forming products of a chains in A connecting 
an index from the set H to an index from the set K 
and passing through the set J. (We remark that [3 may 
be equal to zero. If a = we have a principal minor 
and the discussion does not apply. The term passing 
through the set J as used here includes the possibility 
that no element of the set / appears in a given chain.) 
These products are multiplied by a principal minors 
of Aj and a suitable sign attached. Thus we have 



A R9S = Z±a(hJ l o- 1 ) 



a(haj a o- a )Aj> 



(3.5) 



where en is an element from the set K, An is a minor 
of A, and the sum is taken over permutations of the 
set K and partitions of J. 

With the aid of this result we may generalize lemma 
3.4 to obtain: 

Lemma 3.6: Let AeQ, then the formula (3.5) for A RfS 
has at most one nonzero term. 

PROOF: Suppose for contradiction, that there is 
more than one nonzero term in (3.5). To be specific, 
suppose the factor, a(hiJ\k\)a(hzj2k2) appears in 
one term and the factor a(hiJik 2 )a(h2j2ki) appears in 
another term. Consider the principal submatrix A\L] 
in rows and columns hi, h 2 , k u k 2 , and the q distinct 
elements from the sets Ji, J 2 , Ju an d }%• This matrix 
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is of order p + 4 and contains p + 4 nonzero elements 
above the main diagonal. This contradicts the fact 
that AeQ. 

We may also analyze the cycles of A P A. Let 
Ji, . . ., Jr be increasing multiindices of length p, 
then 



A*(/, 



Jr)=A JlJ2 



± J r Ji 



(3.6) 



is an r-cycle of A P A. It is clear that the interchanges 
of rows and columns used above in passing from A RyS 
to A~ ~ will, when applied throughout the cycle (3.6) 
of A P A, not result in any change in sign, On the other 
hand, it is also clear from (3.5) that the product (3.6) 
consists of cycles of A on the elements of J\ with 
principal minors of A. The principal minors of A are, 
of course, just the one cycles of the several matrices 
A p A. 

Turning now to the case of AeQ we have the ap- 
propriate generalization of lemma 3.5. 

Lemma 3.7: Let AeQ, then every nonzero cycle of 
length ^2 of A P A is a product of 2 -cycles of A and 
principal minors of A. 

Lemmas 3.4 and 3.6 generalize the known formulas 
for Jacobi matrices to all elements of Q. Lemmas 
3.5 and 3.7 apparently give new results even for the 
elements of J>*. The cycles of A P A for AeJ>>i do not 
seem to have been studied. 

4. Semipositive Elements of ^ 2 

The semipositive elements of Q seem to have been 
first studied by Parter in [3]. He was aware of a 
slightly earlier paper of Givens [11], which treats 
Jacobi matrices of semipositive type. But a funda- 
mental study of the Jacobi case contained in chapter 2 
of the beautiful book of Gantmacher and Krein [7] 
seems to be almost unknown. (The only noteworthy 
exception to this statement is in the extensive work 
of Karlin on total positivity. In fact, his forthcoming 
book [8] will contain many results which supersede 
to a considerable extent the work in [7] as well as a 
wealth of new applications.) Our work in this section 
has been influenced primarily by Gantmacher and 
Krein and by Parter. 

Definition 4.1: AeJ?2 is of semipositive type if every 
nonzero 2-cycle of A is positive and a# is real for 
7=1, . . ., n. A is principally positive if A is of semi- 
positive type and ajy^0,y=l, . . ., n. 

In this section we shall restrict our consideration 
to the set Q. The reason for this is simply that the 
author is unaware of any interesting applications of 
the semipositive elements of J>i — Q. We start with: 

THEOREM 4.1: Let AeQ be of semipositive type. 
Then there exists a diagonal matrix D uniquely deter- 
mined by A up to the condition d(D) 4= and the signs 
of the elements such that 

A^D^AD 
is a real symmetric matrix. 



This theorem has an interesting history. Parter 
pointed out in [3] that AeQ of semipositive type is 
symmetrizable, basing his result upon previous work 
by Hearon [6] and Goldberg [4, 5]. Independently the 
present author pointed out in [1] that A is symme- 
trizable and, apparently, first noted that this could be 
accomplished with a diagonal matrix. (In fact, the 
paper [1J gives a more or less explicit formula for D 
when A is real and a similar formula is not difficult to 
obtain for A complex.) It is also possible to show that 
the signs of D can be so chosen that A has nonnegative 
elements except possibly on the principal diagonal. 
This leads one to the following theorem which is a 
transcription for the semipositive case of a theorem 
proved in [1], We do not include a proof. 

THEOREM 4.2: Let AeQ be of semipositive type, 
then: 

ia) for each b ^ min (a jJ5 0) the matrix A + bl has a 
unique positive eigenvalue, p(A -f bl), such that for 
every X in the spectrum, cr(A 4- bl), p(A + bl) > |X|. 

(b) If x is an eigenvector of A -f bl corresponding to 
p(A 4- bl) then there exists a nonsingular diagonal 
matrix D such that x = Dy where y > 0. 

(c) The spectrum of A and the spectrum of every princi- 
pal sub matrix of A is real. 

(d) If any element of A is increased, then p(A) is strictly 
increased. 

In order to explore the spectral properties of ele- 
ments of Q which are not so near the surface, it turns 
out to be useful to single out the class Q"' 1 for special 
study. For this purpose, let us consider the nested 
sequence A°=l, A 1 . . . ., A" = A of principal sub- 
matrices associated with AeQ" A in theorem 2.4. It 
will be convenient to set 



d q =d(A<i), q = 0, 1, . . 
d q (k) = d(A« - kl) , q = 0, 1, 



n, 



(4.1) 



Thus do(k) a 1. 

The theory begins with: 

LEMMA 4.3: //AeQ n_1 is of semipositive type, the 
sequence of polynomials (d (X), d(X), . . ., d n (X)), is a 
Sturm sequence for every real value of X. 

PROOF: Since the polynomial d q (\) is of degree q for 
q=0, . .* . n, we need only show that, if d I) (ko) = Q, 



rfp(Xo)rfp-i(X ) < o. 



(4.2) 



In the present notation the formula 2.7 applied to 
A — XI in place of A reads 

d <iW = (%/g ~ X)d q -i(\) - %_ 1 j p «v 9 _ 1 ^-2(X). 

(4.3) 

If we apply (4.3) for q = p 4- 1. multiply by dp-i(X) and 
set X = X , we obtain 

d/;+i(Xo)dp-i(Xo) = ~ %_ 1 j p %i p _ 1 d|-i0w) ^ °' 

since A is of semipositive type. The lemma is thereby 
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proved unless d p -\{Xo) = 0. Suppose, in fact, this is so. 
Applying (4.3) for q = p and setting X = X () shows 
dp-2(K)) — also. Thus, by the same reasoning, 



dp-dko)- 



:rf 1 (X ) = do(\o) = 0, 



contradicting the fact that doiXo) = 1. It follows that 
dp-i{\o) 4= and the lemma is proved. 

THEOREM 4.4: Let AcQ" -1 be of semipositive type. 
Then the eigenvalues of A q , q=l, . . ., n — 1, strictly 
interlace with those of A q+1 . Hence, in particular, the 
eigenvalues of A are distinct. 

PROOF: By Sturm's theorem the Cauchy index 
I^d Q -i(\)ld Q (\) (see Gantmacher [12], Vol. 2) equals 
V q (a) — V q (b) , where V q (X) is the number of changes 
in sign in the Sturm sequence (d (X) , . . ., d Q (X)) 
and where — oo^a<6^H-°°. In the present case 
we have d p (X) = (— l) p X p -\- . . . for each value of p. 
Thus V q {— oc)=0 and V q (-\-°°)=q. The theorem now 
follows immediately from the meaning of the Cauchy 
index. 

If A$Q is of semipositive type and A^Q n ~\ then 
A may or may not have a multiple eigenvalue. For 
example 

3 -1 -1 



A x 



3 






is of semipositive type and not in Q n ~ 
multiple eigenvalue. But 



and has no 



A 2 = 



f 3 
-1 



1 

2 





has the double eigenvalue X = 2. Thus the fact that 
AeQ"~ l has no multiple eigenvalues if it is of semi- 
positive type does not distinguish this class. We shall 
show now that it is the interlacing of the eigenvalues 
in the sequence A 1 , . . ., A n = A which does dis- 
tinguish the class. 

The following theorem is quite general in its appli- 
cation. 

THEOREM 4.5: The matrix A posses a sequence of 
principal suhmatrices with strictly interlacing real 
roots if and only if there is a sequence of principal 
submatrices whose determinants form a Sturm sequence 
for all real X. 

PROOF: The "if part is the content of theorem 4.4, 
the proof of which did not make explicit use of the 
hypothesis AeQ n ~^. Hence we may confine ourselves 
to the only if statement. 

To be specific suppose J\CJ 2 C . . . Cj n is a se- 
quence of increasing multiindices with J q containing 
q distinct elements from the set J n = S n = (1 . . . n). 



Let Jo be the empty set of indices and assume the 
sequence of polynomials. 



(4, (\), ^(X), . . .,/*, (X)). 



(4.4) 



where A J() (X)=l, is such that Aj p (k) has p distinct 
real roots which strictly interlace with those of Aj p+l (k), 
p=l, . . ., n — \. Denote by X^, i, . . ., X p , p the 
zeroes of Aj (X) in increasing order. For each p we 
require that 

Xp+i, i < Xp, i < Xp+i, 2 < • . . < k PfP < Xp+i, p+i. 

In view of the fact that the leading coefficient of 
A j (X) is (— l) p , we may write 

^W^fi^rX), p = l, . . ., n. 

For fixed r we have 

^ +1 (^rM/ H (Ap,r) 

p+l p-l 

| (Xp, r — X p+U q) |J (X p , r — Xp-i, q) 



(7.4) 

In the first product there are r positive factors and 
p + 1 — r negative factors. In the second product there 
are r— 1 positive factors and p — r negative factors. 
Thus there are 2 (p — r) -h 1 negative factors in the 
product and 

Aj 9+1 (kp,r)Aj p _ 1 (\p 9r )<0 9 

for r=l, . . ., p. It follows that the sequence (4.4) 
is a Sturm sequence. This proves the theorem. 

On the basis of theorem 4.5 we can now investigate 
the consequence of Ae^2 possessing a Sturm se- 
quence of principal submatrices. Let us assume that 
the index j p is adjoined to J p -i to obtain J p . Then for 
AeJ^ 2 the formula (2.3) applied to the sequence (4.4) 
becomes 



A Jp (k) = «Hj p )-k)Aj- jp (\) - £ aUMAj-u^W- 

9=1 

(4.5) 
It follows that 

A Jp (kp-i t r)Aj p _ 2 (k p - ltr ) 

p-\ 

= - 2 ^ph)^J-(jpjq)(^p-X,r)A Jp _ 2 {Xp- U r ). 
q=\ 

Clearly we shall have a Sturm sequence only if each 
of the numbers 



2 a{jpjq)Aj- Upjq) {Xp- U r)Aj p _ 2 {Xp- Ur ) > 0, 



(4.6) 
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r=l, . . ., p — 1. Since the inequality is strict, at 
least one term in (4.6) must be different from zero. 
On the other hand, we have conditions of the form 
(4.6) for the n — \ different values of p, p = 2, . . ., rc, 
and A has at most n — 1 nonzero 2-cycles. The same 
2-cycle cannot appear in the inequality (4.6) for two 
different values of p. It follows that the left side of 
(4.6) contains exactly 1-term for each value of 

p = 2, . . ., n, 

and A has exactly n — 1 nonzero 2-cycles. But this im- 
plies that (4.6) becomes 



a(jpJ P -M) p _ 2 {\ p - Ur ) >0 



(4.7) 



p = 2, . . ., n. Therefore each of the n — \ nonzero 
2-cycles of A is positive. This result is summarized in 
the following theorem: 

THEOREM 4.6: Let ke^ 2 possess a sequence of 
principal submatrices with strictly interlacing roots. 
Then AeQ nl and A is of semipositive type. 

Theorems 4.4 and 4.6 nail down the most important 
spectral property of the semipositive elements of Q n ~ l . 
Since the elements of Q n ~ l are essentially Jacobi mat- 
rices by virtue of theorem 2.6, we can hardly argue 
that theorem 4.4 is a new result. On the other hand, 
theorem 4.6 appears to be a new result and, as it com- 
plements theorem 4.4, we felt it appropriate to include 
both. 

The foregoing analysis can be used conveniently 
to study the spectral multiplicity of the elements Q 
not in Q n ~ l . Our intention is to indicate briefly the re- 
sults of Parter, but by a method somewhat more alge- 
braic in character than his. We only sketch the results 
here, referring the reader to [3] for a more detailed 
argument. 

We use theorem 2.6 to derive a kind of canonical 
form for an element of Q — Q 71 ' 1 . By theorem 2.5 
AeQ — Q nX has at least one row with more than two 
nonzero off-diagonal elements. Suppose the pth row 
is one such with the nonzero elements in columns 
<7i, . . ., q a , a > 2. By a modification of the argument 
used to prove theorem 2.6 there exists a permutation 
matrix P such that the matrix PAP' has the follow- 
ing form: 



PAP' 



/A u 
ri 




C\ 



c 2 



. 


. 


r 2 . 


. r c 
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. 



(4.8) 







C a 







in which An is a square matrix of order m 9 i= 1, . . ., a; 
C\ is a column vector of zeroes except for the /iith 
entry which is a QlP ; r\ is a row vector of zeroes except 
for the nonzero element a p , Ql in the n t th entry; c*, 
i = 2, . . ., a, is a column of zeroes except for the 
nonzero element a q . p in the first position: 



n, i 



= 2, 



tion; and each block An is either a Jacobi matrix or a 
matrix of the form (4.8). Thus P A P' is recursively 
defined in the general case. 

Now we apply the basic recurrence formula (2.3) 
to PA P' —XI to obtain 



d{P A P'-kI) = (a pp -\)d{A u (\) 



d(A aa (k)] 



-£ S(p, qi)d(A H (\))l[ d(Ajj(\)). (4.9) 



is a row of zeroes except for a pq . 4= in the first posi 



j*i 



In (4.9) An is the matrix obtained by deleting the last 
row and column of An if i = 1; otherwise it is the matrix 
obtained by deleting the first row and column of An, 
Since a > 3, it follows immediately from (4.9) that 
AeQ — Q n ~ l has the multiple eigenvalue A if A is an 
eigenvalue of at least three of the submatrices An, 
Suppose, on the other hand, that A is an eigenvalue of 
A of multiplicity at least two. Then A must be an 
eigenvalue of the principal submatrix of A, A(p), i.e., 
Ao must be an eigenvalue of An for some value of i, 
say for i = k. Thus 



= d(PA P'-\ I) = -a(p, q k )d(A kk (\ti)fl ^(Ao)). 

j±k 

Hence, either Ajj(ko) = for some value of j =\= k, or 
(A ichi^o)) = 0. In the latter case the block A k k is not a 
Jacobi matrix and must itself have Ao as a multiple 
eigenvalue. Clearly the argument up to this point can 
be repeated for the block A^k, so we may confine 
ourselves to the case where Ajj(\ ) == for some value 
of j 4= k. Thus after a finite number of steps we must 
arrive at a matrix having the form (4.8) such that each 
diagonal block is a Jacobi matrix. But it is easy to see 
that such a matrix can have a multiple eigenvalue Ao 
only if at least three of the blocks have Ao as an 
eigenvalue. 

Now, if a diagonal block An in the representation 
(4.8) is not a Jacobi matrix, there is some row of A, 
say the ^th row, having at least three nonzero off- 
diagonal elements. This row could have been used in 
the beginning to set up the form (4.8). Moreover, if 
this is done, we see that Ao a multiple eigenvalue of 
A implies A, is an eigenvalue of at least three of the 
diagonal blocks. 

This sketch establishes the theorem of Parter which 
we here formulate quite differently than his original 
graph-theoretic formulation. 

THEOREM 4.7 (Parter): Let AeQ — Q nl be of semi- 
positive type. Then A has the multiple eigenvalue Ao 
if and only if there exists a permutation matrix P such 
that PAP' has the form (478) with a >= 3 and A an 
eigenvalue of at least three of the diagonal blocks A u . 

On the basis of this theorem it is obvious that the 
assertions made regarding the matrices A\ and A 2 in 
subsection 4.2 are valid. 

The theorem on interlacing proved above, theorem 
4.4, for the elements of Q n ~ x is of considerable interest 
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in the applications to mechanics. That it does not 
extend in precisely this form to the elements of 
Q — Q n ~ l is quite clear. Nevertheless one may ask 
when X can be an eigenvalue of AeQ — Q n ~ l and also 
an eigenvalue of a certain type of principal submatrix. 
We shall not pursue this question in the present paper. 
It will be taken up elsewhere when we study the 
application of matrices in ^ 2 of semi-positive type 
to problems in mechanical vibrations. Similarly the 
results on eigenvectors are more strongly motivated 
in a mechanical setting. On the other hand, some of 
the results of section 3 have special application to 
the semipositive elements of J 2 and we take up this 
subject next. 

We assume for the remainder of this section that A 
is positive semidefinite. This will be the case, for 
example, if A is diagonally dominant. Lemmas 3.4 
and 3.5 now imply 

THEOREM 4.8: Let AeQ be of semipositive type and 
positive semidefinite. Then cof A has no zero elements 
and is a matrix of positive cyclic type. In particular, 
if A is positive definite, then A ~ l has no zero elements 
and is of positive cyclic type. 

It should be mentioned at this point that some of 
the theory of matrices of positive cyclic type is worked 
out in [2]. In the literature of mathematical economics 
these matrices are usuaDy called Morishima matrices 
in honor of Professor Michio Morishima who first 
characterized them. The analogue of the Perron- 
Frobenius theorem is valid for such matrices. Since 
the largest eigenvalue of A ~ l is related to the smallest 
eigenvalue of A we have 

Corollary 4.9: Let AeQ be of semipositive type 
and suppose A is positive definite. Then if 



Xi 



k n 



are the eigenvalues of A in increasing order, we have 



0<\i < ki<k n , i = 2, 



n-l 



i.e., Xi and X w are simple eigenvalues of A. Moreover 
if any element of A' 1 is increased, then Xi is strictly 
decreased. 

We feel the relationship between theorem 4.7 and 
corollary 4.9 should be stressed. If AeQ — Q n ~ l , is 
semipositive and positive definite, and if A has a 
multiple eigenvalue X*, then X* cannot be equal to 
X ? or to X„. Thus for the matrix A 2 of subsection 4.2 
the double eigenvalue X = 2 is neither the smallest nor 
the largest eigenvalue. 

The fact that A ~ * is of positive cyclic type enables 
us to say even more about A itself. Every matrix of 
positive cyclic type is similar to a positive matrix via 
a diagonal matrix D. Thus, if the conditions of corol- 
lary 4.9 are satisfied there exists a diagonal matrix D 
such that the matrix B~ x = D~ l A x D satisfies 



B" 1 >0. 



(4.10) 



Now (4.10) implies that B itself is monotone, i.e., if 
Bx > then x > 0. Consequently, we deduce that A 



is similar by the diagonal matrix D to a monotone 
matrix. 

In the special case where A^ ^ 0, i^f and A is 
positive cyclic type, then A is both monotone and of 
positive cyclic type. In this case the inverse of A is 
positive, A _1 > 0. That A _1 ^ follows from the well- 
known connection between monotone matrices and 
nonnegative matrices. In order to conclude that 
A' 1 >0 we must appeal to the cofactor formula (3.3), 
and lemma 3.4. This result is a generalization of a 
theorem of Gautmacher and Krein. 

In conclusion a few remarks on total positivity are 
in order. We cannot use lemma 3.7 directly in the same 
way as lemma 3.5, to conclude that A p A is of positive 
cyclic type, since lemma 3.7 is weaker than 3.5. In 
fact, we have seen above that there exist elements 
AeQ — Q n ~ 1 having multiple eigenvalues. Since totally 
positive matrices enjoy the interlacing property 
characterized in theorem 4.4., such elements of 
Q — Qn-i canno t be totally positive. This question 
therefore remains open. It may be formulated as 
follows: Which elements AeQ of semipositive type 
have the property that A P A is of nonnegative cyclic 
type for p = l, . . . , n?(Oi course, we know that 
every such A must be positive semidefinite.) 



5. Seminegative Elements of J 1 

The seminegative elements of Jt first became of 
interest through the work of Quirk and Ruppert [10] 
in qualitative economics. We restrict ourselves here 
to real matrices in ^ 2 - As usual we also confine our- 
selves to the consideration of elements of Q. 

Definition 5.1: AeQ is of seminegative type if every 
nonzero 2-cycle of A is negative. A is principally nega- 
tive if A is of seminegative type and 

a n ^ °' 7=1. • • •> n. 

The following result is established by an argument 
entirely similar to that used to prove a portion of the- 
orem 2 in [1]. We omit the proof. 

Lemma 5.2: Let AeQ be of seminegative type and set 
A = E + A* where E = diag {an, . . ., a nn } and 
A*={a* k } with af k = a jk if j * k and aJ^O. A* is 
similar to a skew-symmetric matrix by a diagonal 
matrix D, i.e., there exists a diagonal matrix D such 
that 

D^AD^E + B, 

where B is skew symmetric. 

Lemma 5.2 leads immediately to the following re- 
result: 

Lemma 5.3: Let AeQ be of seminegative type and have 
a zero principal diagonal. Then all of the eigenvalues 
of A lie on the imaginary axis in the complex plane. 

We shall conclude the paper by formulating the 
Quirk-Ruppert theorem somewhat more elegantly 
than it appears in [10]; (see, however, [9] where the 
result is also better formulated). 
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Definition 5.2: If A = (a^)" is a real matrix we asso- 
ciate with it the matrix sgn A = (af k )i where a* k = sgn 
ajic, j, & = 1, . . . n, and the equivalence class Qa 
defined by 

Q A = {B\sgnB = sgn/4}. 

The equivalence class Qa may be called a qualitative 
matrix. 

Definition 5.3: A is a sign stable matrix if BeQ A 
implies B is a stable matrix, i.e., all the eigenvalues 
of B lie in the open left half of the complex plane. 

Theorem 5.3: (Quirk-Ruppert) Let A be an irreduc- 
ible matrix. Then A is sign stable if and only if 

(i) AeQ and A is principally negative, 
(ii) ajj < for some value of], and 
(hi) d(A) 4= 0. 
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